Rotationally-invariant slave-bosons for Strongly Correlated Superconductors 
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We extend the rotationally invariant formulation of the slave-boson method to superconducting 
states. This generalization, building on the recent work by Lechermann et al. [Phys. Rev. B 
76, 155102 (2007)], allows to study superconductivity in strongly correlated systems. We apply the 
formalism to a specific case of strongly correlated superconductivity, as that found in a multi-orbital 
Hubbard model for alkali-doped fuUerides, where the superconducting pairing has phonic origin, yet 
it has been shown to be favored by strong correlation owing to the symmetry of the interaction. The 
method allows to treat on the same footing the strong correlation effects and the interorbital interac- 
tions driving superconductivity, and to capture the physics of strongly correlated superconductivity, 
in which the proximity to a Mott transition favors the superconducting phenomenon. 
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I. INTRODUCTION 

The theoretical description of strongly correlated sys- 
tems and of the prototypical models introduced to un- 
derstand their behavior plays a central role in modern 
many-body theory. Even if the number of materials of 
interest in which the mutual interaction between elec- 
trons has been identified as relevant is now countless, 
there is no doubt that the main trigger for the devel- 
opment of the correlated-electron field has been the dis- 
covery of high-temperature superconductivity in doped 
correlated insulators such as the copper oxides. Yet, 
the link between strong correlation and high-temperature 
superconductivity has not been established unambigu- 
ously, which prompts for theoretical methods able to de- 
scribe the superconducting phenomenon in the presence 
of strong electron correlations. 

One of the main reasons why strongly correlated sys- 
tems and their properties are, at the same time, interest- 
ing and hard to solve is that they are intrinsically out of 
weak-coupling regimes, where a perturbative expansion 
can be performed. Starting from an uncorrelated system 
and imagining to continuously increase the degree of cor- 
relation, the relevance of local repulsion gradually intro- 
duces constraints to the electronic motion, leading even- 
tually to the localization of the carriers (Mott transition) . 
Thus a proper method for correlated electrons should be 
able to introduce local constraints onto an otherwise un- 
correlated state, that would be naturally delocalized, i.e., 
spatially unconstrained. 

A popular strategy which formally implements this 
point of view is based on slave bosons^'^. Within these 
approaches, the Hilbert space is enlarged to include, be- 
sides fermionic degrees of freedom associated to Lan- 
dau quasiparticles, suitable extra degrees of freedom 
of bosonic character which are typically related to lo- 
cal states. The auxiliary (slave) degrees of freedom 
are then treated in a mean-field approximation, lead- 
ing to an effective low-energy theory for the quasipar- 
ticles. The high-energy physics can only be recovered 



introducing fluctuations of the fields describing the aux- 
iliary particles^. We can see this strategy as a way to 
enforce a local point of view, which is expected to be 
correct in very strong coupling, starting from delocalized 
non-interacting states. In its most popular version, in- 
troduced by Kotliar and Ruckenstein^., one defines one 
boson for each local configuration (namely, |0), |t), 
ItJ,)), and the equivalence between the physical Hilbert 
space and the new extended space is enforced by impos- 
ing constraints which, as we shall discuss below, imply 
that the local configurations should be coherently labeled 
by the fermionic and bosonic degrees of freedom, and that 
precisely one boson should be present on each lattice site. 

Yet, as thoroughly discussed in Ref. 4, the standard 
4-boson representation^ for the single-band Hubbard 
model, as well as its simplest generalizations to multi- 
orbital models^, are not suitable to handle arbitrary 
forms of the interaction Hamiltonian, characterized by 
terms which cannot be put in the form of density-density 
interactions, such as exchange interactions associated to 
the Hund's rule coupling. Furthermore, even for pure 
density-density interactions, the Kotliar-Ruckenstein ap- 
proach still remains inadequate to handle charge sym- 
metry breaking order parameters such as the supercon- 
ducting one in the Hubbard model. Slave-boson ap- 
proaches to superconductivity in the Hubbard model 
have indeed mostly used the approximately equivalent 
strong-coupling t-J model, and have been based on spe- 
cific assumptions^. 

The first attempt in overcoming the inadequacy of 
Kotliar-Ruckenstein's representation was made by Li et 
alX, who proposed a spin-rotation invariant slave-boson 
formulation of the single band Hubbard model, while 
in Ref. @ Fresard and Wolfle introduced a more gen- 
eral representation for single-band models, in which spin 
and charge degrees of freedom are treated on the same 
footing and rotational invariance involves both spin and 
particle-hole transformations. Although the formalism 
presented by these authors refers only to a 4-state sys- 
tem, with the slave-boson fields labeled in correspon- 
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dence with the specific SU{2) (E) SU{2) generators of spin 
and particle-hole rotations, it already has all the ingredi- 
ents required for describing systems with local supercon- 
ducting pairing. Such method has been indeed applied 
to the single-band attractive Hubbard model in Ref. 9,. 
Developing the ideas of these pioneering works in a more 
systematic way, Lechermann et al.'^ finally built a com- 
pletely basis-independent slave-boson formalism, suitable 
to describe, within a generic multi-orbital model, any 
arbitrary form of local interaction. As it is found for 
Kotliar-Ruckenstein's approach, it is worth mentioning 
that, at mean-field level, such formalism turns out to be 
equivalent^ to analogous extensions of the Gutzwiller 

11112 

approacbiiii^. 

While the possibility of extending the formalism to 
superconducting states is mentioned in Ref. 0, the ex- 
plicit derivation is limited to normal solutions, impos- 
ing no charge symmetry breaking. In the present work, 
instead, we lift this restriction and consider explicitly 
the more general case of full rotational-invariance under 
any local transformation of the electronic degrees of free- 
dom, and we apply the formalism to solve a three-orbital 
model which has been proposed to describe alkali-doped 
fuUeridesi^ii^. Besides its relevance to the fullerides, the 
model has important properties that led us to choose it 
as an optimal benchmark for our method. The model 
has indeed been shown to present "strongly correlated 
superconductivity"—, i.e., the enhancement of phonon 
mediated superconductivity in the proximity of a Mott 
transition. The key of the phenomenon is that a small 
attraction which involves orbital and spin degrees of free- 
dom is not screened when charge fluctuations are frozen 
by strong correlations. This leads to an enhancement of 
superconductivity, since the unscreened attraction now 
acts on strongly renormalized quasiparticles with a larger 
effective density of states. This effect has been identified 
using Dynamical Mean-Field Theory (DMFT)^^, which 
fully takes into account local quantum fluctuations, but 
it has not been reproduced by ordinary slave-boson meth- 
ods due to the difficulties in treating interactions which 
are not of the charge-charge form, such as those driving 
superconductivity in the model we are dealing with. In 
this light, the model is an ideal test ground of the ability 
of the rotationally-invariant slave boson method in ac- 
curately treating general forms of interactions. On the 
other hand, the model has only local (on-site) interac- 
tions, which simplifies the approach. 

The paper is organized as follows. In Sec.[TT]we intro- 
duce the rotationally-invariant slave boson method for 
models with local superconducting pairing. In Sec. IIIII 
we present the multi-orbital model used for the descrip- 
tion of fullerenes, illustrating the way it can be solved by 
means of the slave boson approach. In Sec. II VI we present 
the results obtained with our method, and finally Sec. |V] 
is dedicated to concluding remarks and perspectives. 



II. THE GENERAL FORMALISM 

In this section we will explicitly extend the 
rotationally-invariant slave-boson formalism introduced 
by Lechermann et alA to the possibility of describing su- 
perconducting states. To facilitate the reading and the 
comparison with the formalism of Ref. whenever pos- 
sible we shall use the same notation for corresponding 
quantities. 

A. Motivations 

Without entering in details, we shall first provide a 
brief reminder on slave-boson formulations, in order to 
face the difficulties encountered with non-invariant ap- 
proaches such as Kotliar-Ruckenstein's one. 

In a generic multi-orbital model the local Hilbert space 
of electronic states is defined as the set of all the possible 
"atomic" configurations at a given lattice site (for sim- 
plicity, we will drop site indices throughout this section); 
a natural choice for the basis set of this space is provided 
by the 2'^'' Pock states 

\n) = (4)"' • • • |vac), K = 0, 1] (1) 

where a = 1, . . . , M are the local orbital species and d]^ 
the corresponding electron-creation operators. A slave- 
boson representation is then constructed by mapping the 
local Hilbert space Ti. (e.g., the Fock states |n)) onto an 
"enlarged" Hilbert space "H generated by the tensor prod- 
ucts of boson operators (/ij, and auxiliary fermion opera- 
tors /t; 

n : {|n)} ^ 2£ : {n^ ('^t)'^1vac) ® \n) f] . 

In the above expression \ n) f refers to the Fock states gen- 
erated by the auxiliary fermions , which correspond to 
quasiparticle (QP) degrees of freedom: their presence en- 
sures the possibility of describing Fermi liquid properties 
within the auxiliary-fields representation (note that the 
orbital basis for quasiparticle degrees of freedom may not 
coincide, in general, with that of the physical electron op- 
erator djj). On the other hand, an arbitrary number of 
auxiliary bosons (f)^, for each species /i, can in princi- 
ple be present in the enlarged Hilbert space, unless some 
constraints, which characterize the specific form of the 
slave-boson representation, are imposed to its states. In 
other words, a given representation is defined by the way 
in which the auxiliary states are selected out of the en- 
larged space H in order to represent uniquely the original 
physical states 

\n) I — > \n), 

\n) EE J2 K{n,{N^},m) [] {c^lf'\yac) \m)f. 

(2) 
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Needless to say, the choice of a specific representation 
K, apart from being consistent with the above unique- 
ness assumption, must also provide some simplifications 
in the (local) interaction Hamiltonian: the whole purpose 
of introducing auxiliary bosons is indeed the possibility 
of writing local interactions as a sum of quadratic terms 
in the boson fields, at the expense of a larger number 
of degrees of freedom and a more complex structure of 
hopping terms. 

In multi-orbital generalization of Kotliar-Ruckenstein's 
approach, 2^^ boson fields (j)^ = 0„ are introduced in 
correspondence with the original Fock states |n), and the 
representation of such states in the enlarged Hilbert space 
reads 

\n) = 4|vac) ® \n)f. (3) 

The "physical states" in Ti are therefore those states con- 
taining exactly one boson and whose quasiparticle con- 
tent \n) f matches the Fock configuration associated to 
the boson field 0„. From Eq. ([3]) it should be evident 
the non-invariant nature of this representation under ro- 
tations of the quantization basis. Consider, in fact, an 
SU{M) rotation of the orbital indicesi^ 

4 = T.pUo.0dl, fl ^ Y^pUcpIl- (4) 

This rotation will induce a corresponding unitary trans- 
formation on both physical and QP Fock states, \n) = 
X]m^(C^)nm|w), SO that the representation of physical 
states would now read 

\m) = XI {^lin^'lMnm') |vac) ® |to')/ 

n 'm' 

= Y.^Lm'\^^c)®\^')f- (5) 

ra' 

In the new orbital basis, therefore, the slave-boson repre- 
sentation do not retain its original form, and more specif- 
ically the definite relation between physical states and 
their quasiparticle content no longer holds. As discussed 
in Ref. i8|, indeed, only disentangling physical and quasi- 
particle degrees of freedom it becomes possible to formu- 
late rotationally-invariant slave-boson representations. 

As a consequence of its basis-dependent nature, 
Kotliar-Ruckenstein's approach can be applied only to 
systems whose local Hamiltonian can be written, in an 
appropriate basis, in terms of purely orbital-density op- 
erators fla — dl^da, 

Hloc = e^flg + Waf3 flaflf}, (6) 

a af3 

i.e., when the Fock states \n) are eigenstates of H\oc- 
In this case, the representation of H\oc in the enlarged 
Hilbert space can be easily written as a free-boson Hamil- 
tonian, 

n 



with En = J2a ^a'^a + J2ai3 '^o-P Wc remark, how- 

ever, that the definite relation imposed between quasipar- 
ticle degrees of freedom and the (physical) Fock content 
of boson fields inhibits the development of those spon- 
taneous symmetry-breaking order-parameters that can- 
not be expressed in terms of orbital-density operators 
(e.g., superconductivity, magnetization perpendicular to 
the spin-quantization axis, etc.). 

B. Representation of physical states 

The electron Hamiltonian for a generic multi-orbital 
model with purely local interactions is given by 

H = i?kin+5]i?locM, (8) 

i 

iJkin = XXea/3(k)4a'^k/3, (9) 
k Q/3 

where all the local terms, including the chemical potential 
and the orbital energy levels, are included in -ffioc, so that 
Ekea/3(k)=0. 

In comparison to the previous subsection, we choose 
here, as the basis set for the (physical) local Hilbert 
space, a generic set of states {|^)} (not necessarily Fock 
states) that are eigenstates of the local particle-number 
operator n^''^ = X)q=i d'^a'^a, with eigenvalues Na- Even- 
tually, among these sets, we can choose the eigenstates 
{|r)} of the local Hamiltonian, since iJioc commutes with 
the local number operator. The basis set for quasiparti- 
cle states, instead, is still given by the Fock states \n) f 
generated by the auxiliary fermion operators /j,. 

As discussed previously in pointing out the limita- 
tions of Kotliar-Ruckenstein's approach, the key ingre- 
dient in constructing rotationally-invariant slave-boson 
representations is to disentangle physical and quasipar- 
ticle degrees of freedom'*. Therefore, we introduce a set 
of auxiliary boson fields <j}^ = (pAn associated, in princi- 
ple, to each pair (|A), \n) f) of physical and quasiparticle 
states, without assuming any a priori relation between 
those states in the enlarged Hilbert space representation. 
Depending on the phases one takes into account, how- 
ever, there exist some limitations in the possible 2£-states 
4'\n\^^'^) ®W) f which can figure in the representation of 
a physical state |A), 

lA) « X^^Jvac)® |n)/. (10) 

n 

Indeed, if we limit to normal phases as in Ref. 0, it is 
sufficient to consider, for a given state \A), only those 
states \n) f which have exactly the same number of parti- 
cles of \A); in other words, physical states with a definite 
number of electrons are represented by a superposition 
of auxiliary states characterized by the same number of 
quasiparticles. On the other hand, when allowing for the 
spontaneous breaking of particle-number conservation, as 
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in superconducting states, we need to consider, for each 
state \A), all the Fock states \n) f characterized by 



M 
La=l 



Na 



(mod 2) 0, 



i.e., all the 2^^~^ quasiparticle states with the same 
statistics of \A). While the former representation is in- 
variant only under rotations of the QP basis that are 
block-diagonal in the quasiparticle occupation number 
^^^-^riQ, the latter is invariant under a larger class of 
QP rotations, represented by all the unitary transforma- 
tions that preserve the statistics of quasiparticle states: 
in such representation, the quasiparticle number opera- 
tor is no longer a conserved quantity, and its expectation 
value does not correspond to any physical observable, as 
particle-hole transformations may change its value. 

The explicit representation of \A) will then read, in the 
fully-invariant formalism, 



1 



1 ^ 



'/'Li vac) (E)\n)f, 



(11) 



with the sum running over the 2*^^^ QP Fock states 
whose particle- number parity equals that of \A). It is 
worthwhile to remark that the above representation can- 
not be further enlarged, including, for example, in the 
definition of \A), the remaining quasiparticle states with 
opposite statistics. This, in fact, would lead to unphys- 
ical results such as non-vanishing expectation values of 
odd numbers of fermion operators. 



Constraints 

In order to characterize uniquely the physical states 
among all the states of the enlarged Hilbert space H, it 
is necessary and sufficient that the selected states satisfy, 
as operator identities, the following constraints: 



An 

EE<^L'^A„'(n'|/I/a'ki) 

A nn' 

EE<^!in'/'A.'(n'l/t/I» 



1, 

J aJ a 



(12) 

(13) 
(14) 



The first two types of constraints are already present in 
the normal-phase formalism of Lcchcrmann et al.^: the 
first equation, indeed, limits the physical subspace of ?£ 
to one-boson states only, while the second set of con- 
straints ensures the rotational invariance of (jlip under 
QP rotations that preserve quasiparticle number. On 
the other hand, the last set of constraints promotes the 
rotational invariance to particle-hole rotations, enabling 
the non-conservation of quasiparticle number. 

It is worthwhile to remark that, for single-band models 
(M = 2), the above equations reduce to the same set of 



constraints characterizing the SU{2)(E)SU{2) spin-charge 
invariant formalism introduced in Refs. 18| andl9j, as long 
as the appropriate changes in notation are made. For 
this purpose, Ref. 9 provides a useful link between our 
notation and the one presented in Ref. '3. 

A more compact form of Eqs. p2j|14p will be derived 
in Sec. IIIDl where we shall relate the role of constraints 
to the gauge group structure of the slave-boson represen- 
tation. 



C. Physical electron operator 

The representation of the physical electron creation op- 
erator in the enlarged Hilbert space is defined by 



4l^)=^(A|4|i?) \A). 



(15) 



When the constraints (|12m4p are satisfied exactly, its 
expression in terms of bosons and quasiparticle operators 
reads 



1 

M 



AB, nm, 13 



Brn X 



X [{n\fl\m)fl + {n\Mm)f0], (16) 

with the normalization factor 1/M coming from the fol- 
lowing relation: 

E [i^\fl\p)fl\p) + in\fp\p)U\p)] = M\n). (17) 

We can thus summarize the non-diagonal relation be- 
tween physical and quasiparticle degrees of freedom in 
the form 



£ = Ri'iw4 + Ri'Mf0 (18) 



(summation over repeated indices is implied), where we 
have defined the i?-matrix operators as 



R 



ip) 



^An 



AB . 71771 



1 

M 



E 

AB. nrn 



(19) 
(20) 



In the above expressions, we have taken advantage of 
the reality of the matrix elements between Fock states, 
(n|/^|m) = (ml/^ln), in order to guarantee the correct 
transformation properties of the i?-operators under the 
gauge group transformations discussed in Sec. IIIDl 

At the saddle-point level, on the other hand, when the 
boson fields are treated as probability amplitudes and the 
constraints are satisfied only on average, the expression 
of dj, must be modified^, in order to recover the cor- 
rect normalization of transition amplitudes in the non- 
interacting limit. For this purpose, it is easier to define 
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the physical electron operator in the orbital basis {A} in 
which the quasiparticle and quasihole density matrices 
are diagonal, 



Ap) 



(t>*An(l}Arn(m\flff3\n) 

Anrn 
A 

Anrn 



(21) 



(22) 



where ^\ = X^An I^AnP'^A is the probability to find the 
system in a state such that nx = I, i.e., with a quasi- 



particle in the orbital A (note that J^An 



1). In 



these expressions, the quasiparticle operators referred to 
the new orbitals are related to the old ones by the unitary 
transformation 



(23) 



while the boson fields transform with the corresponding 
rotation of the Fock states (summation over repeated in- 
dices is implied) 

\n)f=U{U)„m\m)^, (j)An=U{U)nm^Am- (24) 

In such basis, the transition amplitude between states 
with r^;^ = in the initial [final] configuration, and n\ — 1 
in the final [initial] one, must be normalized by the factor 
1 / •\/^a(1 — yielding the following expression for the 
physical electron operator: 



^1 - E 



7n 



X 



AB.n.n,xV^m^^Im) 

X [^(nlTAllm)-!/)]^ -I- {n\iljx\m)'4!\ 



(25) 



Rotating back to the original basis, we finally get, for the 
saddle-point expressions of the i?-matrices, 



R' 



a(3 y 



AB, nm, 7 

E (^Mii^) 

AB, nm, 7 



where 



i„(l>Brn{m\f'rHMjf}, 

in<l>Br,i{m\fl\n)Mp^, 
(26) 



(27) 



7/3 



is the particle-hole symmetrized version of the normal- 
ization factor, expressed in the original basis. 



D. Functional integral representation 

The partition function of a generic multi-orbital Hamil- 
tonian ([5]) can be formally written, in terms of auxiliary 
fields, as 



2 = Jv[fJ^]vm,{A}] 



(28) 



where we have introduced, along with slave bosons and 
auxiliary fermions, a set of Lagrange multiplier fields 
{^i(r)} that allow to enforce, at each lattice site i and 
imaginary-time value r, the constraints p2j|14p . The La- 
grangian functional entering the above expression reads 



H 



loci 



H 



kin 



(29) 



(except for lattice sites, summation over repeated indices 
is always implied throughout this section), where H^^^ 
and are, respectively, the representatives of the lo- 
cal and kinetic part of the Hamiltonian ([8]) in the enlarged 
Hilbert space, while i?const contains the constraint inter- 
actions between auxiliary fields and Lagrange multipliers. 

In order to derive the expressions of the Hamiltonian 
terms in (j29p . and thereby identify the underlying sym- 
metry group of the Lagrangian, it is convenient to collect 
all the local fermionic degrees of freedom (either physi- 
cal or auxiliary) into a 2Af -component Nambu-Gor'kov 
spinor: 



{da,i} 



In such formalism, the representation of physical elec- 
trons in terms of bosons and quasiparticles (fT8|) is simply 
written as = Ri^i, where the local 2M x 2M matrix 
operator 



R» = R[0,] 



R(p)[0,] R('')[0,]* 
R('')[0,] R(P)[0,]* 



(30) 



is defined in terms of the boson fields (t>An,i associated to 
the corresponding site i (note that, to lighten the nota- 
tion, site indices were omitted in previous sections). The 
representation of the kinetic Hamiltonian is then readily 
obtained as 



ij 

i^*|R|t,,R,vi/^, 



ty 

-t* 



(31) 
(32) 
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tij being the real-space hopping matrix, whose Fourier 
transform gives the band dispersion matrix e(k) defined 
in Eq. ^. 

On the other hand, the local terms of the model Hamil- 
tonian, which may include any kind of on-site interac- 
tion between (physical) electrons, are represented, within 
the enlarged Hilbert space, by a purely quadratic boson 
Hamiltonian: 

^i„Jz] = {A\H,,,\B)^\^^^^Bn. = Er4,,^^<f>Tn,, (33) 
where the physical states {|r)} denote the eigenstates of 



Gauge invariance 

In order to discuss the symmetry structure of the La- 
grangian ([M]) . we begin to notice that the auxiliary- fields 
Hamiltonian — iikin + Ei S.\oc [*] invariant under the 
following gauge transformations: 

*,(t) -> U,(t)*,(t), (34) 
<l>AnAr) ^ e'^"^^^ xUiVU'cl^An'Ar)- (35) 

The unitary matrix Ui (r) in p4p represents an arbitrary 
SO{2M) rotation of quasiparticle operators acting inde- 
pendently on each site, and it is conveniently parameter- 
ized as 



(36) 



the T'^ matrices being a 2 Af -dimensional representation 
of the M{2M — 1) group generators. We note, however, 
that such matrices are not expressed in the usual form of 
an orthogonal group generator (namely, as purely imagi- 
nary antisymmetric matrices), but are instead of the form 



-{Tir -(T|) 



(37) 



with 



and 



denoting, respectively, M x M Her- 
mitian and antisymmetric matrices. In other words, the 
2M-dimensional Nambu spinors do not transform with 
real orthogonal matrices under SO{2M) rotations, even 
though they clearly form a real representation of the 
gauge group: 



U* 



1 

1 



(38) 
(39) 

(40) 



The transformation law of the boson fields, on the other 
hand, is characterized by the unitary transformation 
of Fock states \n) — > Z^[U]„„'|n') that is associated to 
the SO{2M) rotation of quasiparticle operators, plus an 
additional U{1) factor e'^* f"^' under which the Nambu 



spinors are neutral. Following the exponential parame- 
terization of \Ji (t) , we can then similarly write 



(41) 



where jTjJ'jj/ are the SO{2M) group generators expressed 
in the Fock space representation. 

While the gauge invariance of H.^^^ [i] follows immedi- 
ately from its definition ((33)) . in the case of -ffj^jn we still 
need to establish the transformation properties of the 
Ri operators, defined in Eq. ([50)1 . These, however, can 
be readily obtained by writing such operators directly in 
terms of the physical and auxiliary Nambu spinors, in 
matrix notation: 



R,(t) ^ R.(r)U|(T). 



(42) 
(43) 



The above transformation law ensures the gauge invari- 
ance of the physical electron operator Sj = R^^'i, and 
hence of the whole kinetic Hamiltonian ([51]) . 

In the discussion of the gauge invariance of the La- 
grangian (|29p . we are thus left with the time-derivatives 
and constraints terms, whose transformation properties 
are closely related to each other. The time-derivative 
terms, in fact, are clearly not invariant under the trans- 
formations ([34]) and p5|) , which generate inhomogeneous 
terms proportional to the time derivatives of the rota- 
tion parameters (e.g., dr^f and dr^f). Such terms, how- 
ever, can be reabsorbed by a corresponding inhomoge- 
neous transformation of the Lagrange multiplier fields^iii 
(which may be regarded as "gauge bosons" ) , making the 
whole Lagrangian gauge invariant. 

To show how this mechanism works, we rewrite the 
A'P + M{M - I) local constraints US]) and HI]) in the 
following way: 



(44) 



where we have made use of the "orthogonality" between 
the S0{2M) T° matrices. It is then straightforward 
(though somewhat lengthy) to verify that the Fock-space 
generators , introduced in Eq. (j4T|) , are represented 
by 



(45) 



(in other words, the above matrices provide a faithful 
representation of the S0{2M) Lie algebra), so that we 
can finally write: 



^An.i 



1 



^ nn 



(46) 



Together with the time-derivative terms, the above inter- 
actions may be arranged in "covariant derivatives" acting 
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on the auxiliary fields (fermions and bosons), 



(47) 



'An,i 



(48) 



where the role of gauge fields is played by the Lagrange 
multipliers Af and A*^ . It is then easily checked that the 
Lagrangian 



Kloc 



-A" 



(49) 



is indeed invariant under the SO(2M) (g) U{1) gauge 
group, provided the Lagrange multipliers transform as 
gauge fields in the adjoint representatioi*i^: 



AfT" 



A^-idrC^r), (50) 
U,(r)KT"+5,]U|(r). (51) 



Note that the transformation law of ^"T"^ induces a cor- 
responding transformation of A^J^^,^, that has exactly 
the same structure of ([STjl. with the Ui matrix replaced 
hyU[lJ]nn' ■ More precisely, both transformation laws de- 
scend from that of the Lagrange multiplier field ^f(r), 
which for infinitesimal rotations transforms as 

At^At + r'^AUt - ^^r^: + Oie). (52) 

jrabc ]-,gijjg structure constants of the group. 

Finally, it is worth mentioning that for M — 2 (single- 
band models) the orthogonal group 5*0(4) is locally iso- 
morphic to SU (2)(g)S'[/(2), so that the gauge group struc- 
ture of the present formalism reduces to that of the spin- 
charge invariant formalism of Ref. H. 



Gauge fixing 

As discussed in Ref. 8, the gauge invariance of the 
functional integral representation causes Eq. (f28|) to con- 
tain integration over spurious degrees of freedom, namely 
the physically equivalent field configurations connected 
to each other by gauge group trajectories. It is thus nec- 
essary to impose a "gauge fixing" condition that removes 
the integration over the unphysical degrees of freedom, 
as it is usually done in gauge field theoriesi^i. 

We choose to work in the so-called "radial gauge"—, in 
which the complex boson fields are represented through 
real amplitudes and complex phase fields. In this rep- 
resentation, the SO{2M) (g) U{1) gauge transformations 
allow to remove M{2M — 1) -f 1 phase variables, so that 
a corresponding number of boson fields can be reduced 
to purely real amplitudes, with no phase fluctuations. 
The boson fields that remain fully complex, on the other 
hand, continue to display some phase dynamics, which is 
responsible for the incoherent features of the spectrum^ ^ 
(e.g., lower and upper Hubbard bands). 



It is beyond our purpose to enter into the formal details 
of the radial gauge representation, thoroughly derived in 
Ref. [l3- Nevertheless, it is worthwhile to observe, here, 
that such gauge fixing procedure allows to avoid Elitzur's 
theorem^-, which would prevent the slave bosons to ac- 
quire a non-zero expectation value, making therefore le- 
gitimate the use of the saddle-point approximation. 



E. Saddle-point approximation 

The saddle-point approximation of the functional inte- 
gral (p8| is obtained by considering the slave bosons and 
Lagrange multipliers as static variables, corresponding to 
the time-averages of these fields. Moreover, we will as- 
sume a homogeneous spatial structure of the saddle-point 
solution, so that we can finally set: (j)An,i{T) i~+ ^fAm 
A°:{t) ^ A" and A°{t) ^A°. 

In such approximation, all the bosonic amplitudes (pAn 
are assumed to have a constant phase, in contrast to the 
radial gauge representation, where we are allowed to re- 
move (fix) a limited number of complex phases (namely 
M(2M - 1) -I- 1). The phase fluctuations of those flelds 
that remain intrinsically complex are thus ruled out, 
precluding the possibility of describing the high-energy 
physics of a given model. The low-energy features, on 
the other hand, will be suitably described in terms of 
coherent Landau-Bogoliubov quasiparticles, providing a 
Fermi-liquid description of metallic and superconducting 
states. 

At the saddle point, the free energy per site = 
— {1/ PN)hiZ is obtained as the (minimum) stationary 
value of the following free-energy functional: 



n[{cp},{A}] 



\nZf-A° 



(53) 



f3N 

+ J2 V*An[iA\Hloc\B)S^^ 
AB nm 

+ {A°Snm + A^'Jnm) ^ Ab\ ^Bm, 

where N is the total number of sites and Zf represents 
the Gaussian integral over the auxiliary fermions. 



J exp -^J^^dT J2 *lc [dr + h(k)] *k 



(54) 



In the last expression, h(k) denotes the momentum-space 
quasiparticle energy matrix: 

h(k)=RtM(^<Q^) _^^\^^^R[^]+A^T\ (55) 

To evaluate the functional integral in Eq. ([5^ . we note 
that, up to irrelevant constants, Z| = det [dr + h(k)], so 
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that we can write 



InZf = -lndet[a^ + h(k)] 



1 + e-^^'O') 



(56) 



The saddle-point equations are then obtained by set- 
ting to zero all the partial derivatives of the free-energy 
functional (j53p with respect to the slave boson ampli- 
tudes and Lagrange multipliers. For practical calcula- 
tions, however, it is useful to consider a different basis 
set for the M{2M — 1) Lagrange multipliers A"", belong- 
ing the adjoint representation of SO{2M). In place of 
them, in fact, we can equivalently use the set of indepen- 
dent matrix elements of A'^T"-, parameterized as follows: 



Q/3 
Q/3 



A 



A 

n* 



n 

-A* 



(57) 



With this choice, the saddle-point equations can be de- 
rived without knowing the explicit expressions of the 
T" matrices, differentiating the free-energy functional di- 
rectly with respect to Aap and Tlap- To this end, we note 
that ri[{iy9}, {-4}] can be easily expressed in terms of the 
new set of Lagrange multipliers by means of the following 
relation: 



A a /ja 
'Jnn' 



i(n'|1'^(yt"T'^)^'|7i) 



Kp{n'\{fif(,^ fpfi)\n) (58) 



na/3(n'|/t4|n)+n;^(n'|//3/a|n> 



F. Green's functions and observables 

After solving the saddle-point equations, we can finally 
obtain the expressions for the quasiparticle and phys- 
ical electron propagators, conveniently written here in 
Nambu notation. For quasiparticles, the propagator is 
defined as 



D/(k,r) 



(59) 



G/(k,r) 
F/(k,r) 



FUk,-r) 



-GJ(- 



k. 



where G/(k,r)„^ = - {T h^ir) flp{Q)) and i^/(k, t)„^ = 

~ {T f'L\^a^T) fl_i^(Q)) are the normal and anomalous quasi- 
particle Green's functions. Following Eq. (|54p . we can 
then readily write the quasiparticle inverse propagator 
as 



D7i(k,c^) 



h(k). 



(60) 



The physical electron propagator, on the other hand, is 
defined by 



D,(k,T) = -(TSk(r)4(0)), 



(61) 



where ^j^. — (^{^kQ,}, {d-ka}^ is the Nambu spinor con- 
taining the physical degrees of freedom, represented in 
terms of slave boson amplitudes and quasiparticles by 
= R[(^]4'k. The expression for the inverse physical 
propagator is thus written as 



D-i(k,c.) = [Rt]-i[c.-h(k)] R-i. 



(62) 



Using the corresponding expression for the "bare" phys- 
ical propagator. 



Ddo(k,w) = ^-[eo + e(k)]. 



eo -h e(k) = 



Co + e(k) 








-[eo + e(-k)]' 



(63) 
(64) 



(eo represents the one-body part of -ffioc), we can then 
find the saddle-point approximation for the self-energy: 

= c^(l-[RRt]-i)-go 



fR 



ti-i 



n 



n* 



, R-^ (65) 



The cancellation of the fc-dependence, in the above equa- 
tion, follows from the definition of the QP energy matrix, 
Eq. ([55]) . Indeed, this form of the self-energy is just the 
one we would expect from the saddle-point (mean-field) 
approximation, which freezes spatial and dynamical fiuc- 
tuations. From the linear term in w, one readily obtains 
the matrix of quasiparticle spectral weights: 



5Sd 

duj 



UJ=0 



RRK 



(66) 



We conclude this formal section by writing the rep- 
resentation, in the enlarged Hilbert space, of the local 
physical density operator, whose expectation value de- 
fines the average number of electrons per site: 



(67) 



As mentioned previously, we remark that this expression 
differs substantially from that for the local quasiparticle 
density. 



E -^a-^" = E E '^L'/'An, 



(68) 



An a 



which is not a physical quantity and depends from the 
choice of the QP basis set. More generally, the represen- 
tation of any (local) two-particle physical operator may 
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be easily obtained in terms of boson operators only: for 
particle-hole operators we find 

= ^(A|4d^|i3)^0L<^Bn, (69) 

AB n 

while for particle-particle operators we have 

£4 = Y.{A\did\\B) ^ 4,>Bn. (70) 

AB n 



III. APPLICATION TO THE THREE-ORBITAL 
MODEL FOR FULLERIDES 

A. The model 

In this section we present an explicit application of the 
rotationally invariant slave boson approach. To this aim 
we have chosen a three-orbital Hubbard-like model that 
has been used to understand the role of strong correla- 
tions in alkali-doped fuUerides. The physics of alkali- 
doped fuUerene systems AnCgo represents an optimal 
playground in understanding the key role of strong cor- 
relations on high-temperature superconductors: indeed, 
these systems display a relatively high Tc compared to 
ordinary Bardeen-Cooper-Schrieffer (BCS) superconduc- 
tors, and, similarly to what is found in cuprate systems, 
the enhancement of Tc seems to be closely related to the 
proximity of a Mott insulating phas o^^i^^i^"^ . Although 
the nature of the pairing mechanism is most likely differ- 
ent from that characterizing cuprate superconductors, in 
fuUerene systems being due to ordinary electron-phonon 
(vibron) interactional^ and in spite of the different sym- 
metry of the order parameter (s-wave in fuUerides), the 
phase diagram as a function of the inter-molecule separa- 
tion in these systems presents strong similarities to that 
of cuprates as a function of doping^'^, providing an inde- 
pendent example of the key role of electronic correlations 
in enhancing superconductivity. 

The local Hamiltonian describing the Ceo"" molecular 
ion, assuming rotational invariancc within the threefold 
degenerate level <i„ hosting the valence electrons pro- 
vided by the alkali metals, can be written, for a generic 
site i, as 



U.2 r 



2 Si • S,; + — Lj • lii 



(71) 



where the three terms represent, respectively, the global 
on-site Coulomb repulsion, Hund's rules splitting ( Jh > 
0) and the Jahn- Teller coupling between electrons and 
the vibrational modes (vibrons) of Cgo- In the above ex- 
pression, hi = X^ao- aa'^i-o.cr IS the local electron number 
operator (a = 1,2,3 labels the tiu orbitals and a =t,i 
the spin components), while and are, respectively. 



the local spin and orbital angular-momentum operators, 

= ^ ^ dl^^&aa'di^aa', (72) 



Li — ^ dl^^^iabdi^ba, 



(73) 



ab, a 



where (To-o-' are the Pauli matrices and (!>^^ — isbac the 
0(3) group generators characterizing orbital rotations. 

The Jahn- Teller Hamiltonian involves both electron 
and vibron field operators, but if we are interested only in 
the electron dynamics we can formally integrate out the 
vibronic degrees of freedom, obtaining an effective action 
for the electron operators only. If performed exactly, this 
procedure would clearly generate non- instantaneous (i.e., 
time-dependent) interaction terms, preventing a purely 
electronic Hamiltonian formulation of the effective ac- 
tion; however, if we assume the vibronic frequencies to 
be much higher than the relevant electronic scales, we 
can take the anti-adiabatic limit of the electron effec- 
tive action, neglecting retardation effects and consider- 
ing an instantaneous interaction term which preserves 
the symmetries of the original local Hamiltonian. While 
this approximation may be questionable for fuUerides, it 
will not affect our claims since the neglect of retarda- 
tion can only disfavor superconductivity, and analogous 
results have been obtained in a similar model that takes 
into account the phonon dynamics^^. The Jahn- Teller 
interaction can then be reabsorbed in Hund's term, and 
the resulting Hamiltonian is simply given by the first two 
terms of Eq. (|7ip with a renormalized Hund's coupling 
Jh i-^ — J = Jh — |£^jt < 0, i?jT being the char- 
acteristic Jahn- Teller energy gain^^. The net effect of 
the electron- vibron coupling is therefore that of revers- 
ing Hund's rules, favoring atomic configurations with low 
spin and orbital angular momentum. The inversion of the 
Hund's rule is experimentally confirmed by the low-spin 
state of both tetravalen1j2l and trivalent^ fuUerides and 
by the presence of a spin gap^^. 

Given the local Hamiltonian for the Ceo"^ molecular 
ion, the expression for a tight-binding electronic Hamil- 
tonian describing AnCgo solids will then read 

iJlatt = E AaAjM + E W - M E (74) 

ij, ab, a i i 

where i^j* are the inter-site hopping amplitudes (includ- 
ing possible inter-band hybridization terms) and /i is the 
chemical potential controlling the average electron den- 
sity. We should note, however, that inter-band hybridiza- 
tion can actually be avoided whenever the hopping terms 
are restricted to nearest-neighbors only, t^j* = i°^(5^y^ ; in 
such case, indeed, we can exploit the 0(3) orbital sym- 
metry of H\oc in order to diagonalize t"'', so that we can 
set, without loss of generality, t^j" = t°-5ab5(ij)- Through- 
out our analysis we will use the latter expression for the 
hopping matrix elements, focusing in particular on the 
orbitally degenerate case = —t (we will consider the 



10 



possibility of a crystal-field splitting of the three bands 
and of different bandwidths in a future study). 



B. Slave-boson representation of the model 

As pointed out in previous studie a^^i^^'^^i"^" , the zero- 
temperature phase diagram of the model ifTl]) as a func- 
tion of the ratios J/W and U/W (W being the non- 
interacting bandwidth) displays several interesting fea- 
tures, the most striking one being undoubtedly the pres- 
ence of a strongly enhanced superconducting phase in the 
proximity of the metal-to-insulator Mott transition. The 
model represents therefore a valuable test for the rota- 
tionally invariant slave boson method, at the same time 
providing an analytical tool to treat strongly correlated 
superconductors. 



Slave-boson amplitudes 

From the rotationally-invariant slave-boson represen- 
tation of physical states defined by Eq. pT|) . it should 
be clear that there are, in principle, 2^^~^ (with M = 6 
in our model) slave-boson fields (j)An describing the sys- 
tem. However, we must note that when the partition 
function is approximated at saddle-point level, i.e., when 
the slave-boson fields are replaced by their mean-field ex- 
pectation values {(f) An {4>An) = ^An), thc numbcr of 
independent slave-boson amplitudes entering the saddle- 
point equations becomes much smaller, as we will show 
in the following, its specific value depending on which 
symmetries of the model Hamiltonian remain unbroken. 

As discussed at end of Sec. lIII A] we choose the hopping 
matrix to be diagonal and degenerate in both spin and or- 
bital indices, so that the atomic SU{2)®0{i) symmetry 
characterizing H\oc can be promoted to a global symme- 
try for the full lattice Hamiltonian i?iatt- If we impose 
this symmetry to be preserved at saddle-point level, i.e., 
we do not allow for any spin and orbital ordering, we 
must then set to zero all the possible order parameters 
which are not invariant under SU{2) (g) 0(3), and this 
will strongly limit the possibility of independent slave- 
boson amplitudes. Indeed, considering the normal state 
solution of the saddle-point equations, i.e., do not al- 
lowing for a superconducting order parameter, it is quite 
straightforward to prove, using Wigner-Eckart's theorem, 
that the non-zero slave-boson amplitudes must be of the 
form 



<Prn 



{n\T) ^Er), 



(75) 



where we have taken as the basis set for the local physical 
states the eigenstates {|r)} of iJioc, with eigenvalues Ep- 
Note that in this case the quasi-particle Fock states \n) 
have exactly the same number of particles of the physical 
state |r) to which they are linked, assuring the solution 
to represent a normal state; indeed, as long as the latter 



condition is satisfied, no superconducting order parame- 
ter can be ever developed, as can be easily seen taking 
the expectation values of Eq. ([70]) , 



.4) 



(76) 



rr' 



On the other hand, if we do allow for a superconducting 
symmetry-breaking, we must add to the normal ampli- 
tudes defined in Eq. ([75)1 also those amplitudes connect- 
ing physical states to QP states with a different num- 
ber of particles, so that particle number would no longer 
be conserved; however, if we still want to preserve the 
SU{2)(^0{3) symmetry as in the normal case, we should 
consider only those amplitudes which correspond to an 
invariant, with respect to spin and orbital rotations, su- 
perconducting order parameter. Assuming pairing to be 
purely local, corresponding to an s-wave order parame- 
ter, the only invariant pairing amplitude is then given by 
the spin and orbital singlet channel 



t.aj t. 



(77) 



At this point, it is worthwhile to observe that ((77)) repre- 
sents the most favorable pairing channel even if we do not 
explicitly impose the SU (2) (g) 0(3) symmetry, since the 
local pairing attraction is driven by the reversed Hund's 
term, which favors the formation of two-particle states 
locked in the L = S = spin-orbital configuration. Our 
assumption of preserving the SU{2) (g) 0(3) symmetry is 
then fully justified whenever the system turns out to be 
superconducting, since any rotational symmetry break- 
ing pairing would be ruled out by the singlet channel; 
on the other hand, in our study we will only compare 
the superconducting solution with a rotationally invari- 
ant normal state, and therefore we cannot exclude the 
possibility that some other ordered phase would win the 
competition for the lowest-energy phase. 

We can now turn to the problem of establishing the in- 
dependent slave-boson amplitudes required by our model 
in order to describe a superconducting solution charac- 
terized by the L = S = order parameter defined in 
Eq. ([77)) . Using the Wigner-Eckart theorem as for the 
normal state solution, in this case we find the following 
expression for the non-zero amplitudes: 



E 

9=1 



(n|v«|r) 



T|(^t)9^g|r) 



--^iEr,2q) 



^iEr,~2q) 



ipvn = {n\T)(^{ET 
, (n|(^t)9|r) 

T|v«(V't)ir) 

(78) 

We have denoted with $(i?r) the "normal" slave-boson 
amplitudes, which relate physical and quasi-particle 
states characterized by the same number of particles, 
and with ^'(-Er, AiV) the "anomalous" ones, in which the 
number of particles characterizing the QP state \n) differs 
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by AA'' from that of the physical state |r). We remark, 
however, that the presence of non-vanishing anomalous 
amplitudes is not sufficient, by itself, to assure a super- 
conducting solution, the latter requiring the supercon- 
ducting order parameter V'scj which is a specific quadratic 
form in the slave-boson amplitudes, to be finite. The 
normalization factors for the anomalous amplitudes, in 
Eq. (j78p . are chosen in order to simplify the expression 
for the probability associated with the physical configu- 
ration |r), 

n 

= mEr)\' + J2\'^iEr,AN)\\ (79) 

Using Eq. ([75]) . which relates all the slave-boson am- 
plitudes iprn to the independent variables ^(Er) and 
^(E'rj A7V), we are almost ready to write the free-energy 
functional ([55)) in terms of only <i>'s and ^''s, obtaining 
therefore a much smaller number of saddle-point equa- 
tions to be solved. The last step needed to achieve this 
goal, in fact, is to evaluate the matrix elements (njOrflr) 
between the eigenstates of -ffioc , which form the basis for 
the local physical Hilbert space, and the Fock states \n)d 
expressed in terms of the physical electron operators djj^ . 

As discussed in Sec. IIII Al the effective local Hamilto- 
nian for the electron dynamics of the Ceo"" ion is given, 
in the anti-adiabatic limit, by 



U 



loc 



(n-3)2 



fin 



(80) 



2S-S + -L-L 



6(- 



3)^ 



where we have included also the chemical potential term, 
as required by the general formalism described in Sec. [ill 
and we have dropped, for simplicity, the redundant lat- 
tice site index. In comparison to Eq. (|7ip . the Coulomb 
interaction is here written in a particle-hole symmetric 
form by properly redefining the chemical potential 

We can then readily identify the eigenstates of Hioc 
among the atomic multiplets |r) which are simultaneous 
eigenstates of the density operator n and of the orbital 
and spin angular momentum operators and S^, with 
eigenvalues 



Er = E{n,e,s) 

= - 3)^ - 



(81) 



-I- J [2s{s + 1) + i£(^ + 1) + |(n - 3)2] . 
The degeneracy associated to each eigenvalue is given by 

= (2^+l)(2s + l) (82) 

and it is therefore natural to choose as a basis set for 
the corresponding degenerate subspace the simultaneous 



eigenstates of both one of the components of L and S, 
say and 5*2, so that we can finally set 



|r) = |n, (^,4),(s,s.)) 



(83) 



These states must then be expressed in terms of the 
Fock states \n)d^ in order to evaluate the matrix elements 
which characterize Eq. ([78| . We note, however, that 



|n),^n n (4.)"""|vac) K.=0,1] (84) 

a=l a=-\,i 

are not eigenstates of any of the orbital angular mo- 
mentum operators L^, making the representation of the 
atomic multiplets ([83[) in terms of such states a bit in- 
volved: it is more convenient, instead, using the rota- 
tional invariance of H\oc, to choose an orbital basis for 
the physical electron operators in which is diagonal. 



''mcr ~ / ^ '-'ma'J-aa, 



(85) 

[to = 1,0,-1] (86) 



so that the corresponding Fock states 

m, (T 

= U{U)nn'\n')d 



vac) 



(87) 



are eigenstates of both Lz and Sz. Since we are assuming 
the three bands to be degenerate, the rotation of the 
orbital basis ([85[) does not change the form of the kinetic 
term in -ffiatt, while the expression for the singlet pair- 
creation operator tp = c?I|^lx reads, in the new basis. 



'-iT''-U 



(88) 



The representation of the atomic eigenstates jP) in 
terms of the Fock states \n)c is listed in Table [H where 
they are classified according to the quantum numbers 
(n,£,s) which determine, through Eq. ()8ip . the corre- 
sponding eigenvalues; note, however, that for a given 
value of the particle number n, the Pauli-principle pre- 
vents the orbital and spin angular momenta £ and s 
to take independent values, so that each degenerate- 
multiplet can actually be identified specifying only two 
quantum numbers, n and £. For each multiplet, we have 
written out explicitly only the component characterized 
by the maximum value of £z and Sz , all the other compo- 
nents being easily obtainable from the former by repeat- 
edly acting on it with the lowering operators 

L- = (^clc,^ + c^,c„,) , (89) 

^- = E^^i^^T- (90) 
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|r> 


n {I, (s, Sj) 




0,0,0) 


(0, 0) (0, 0) 


$(0,0) 
*(0,0;2) 
*(0,0:4) 
*(0, 0; 6) 


T,0,0> 


1 (1, 1) (i 5) 


$(1,1) 
*(1, 1; 2) 
*(1, 1;4) 


in, 0,0) 


2 (2, 2) (0, 0) 


$(2, 2) 
*(2, 2: 2) 


|T:T,0> 


2 (1: 1) (1, 1) 


$(2, 1) 

^ (2 1" 2) 


4. hT,o,i) - |i,o,T) - |o,Ti,o)l 

V L J 


2 (0, 0) (0, 0) 


$(2,0) 
*(2, 0: -2) 

*{2, 0;4) 


in,T.o) 


3 (2, 2) (i, i) 


<^(3,2) 


^ [lT,Ti,o) + |Ti:0,T>] 


3 (1. 1) (I, i) 


*{3, 1; -2) 
^'(3 1' 2) 


IT,T,T> 


3 (0,0) (|, i) 


$(3,0) 


in,n,o) 


4 (2, 2) (0, 0) 


$(4, 2) 
*(4, 2: -2) 


in,T,T) 


4 (1, 1) (1, 1) 


$(4, 1) 
'I'(4, 1; -2) 


[li-Ti,T> - IT.na) - ITi,0,Ti>] 


4 (0, 0) (0, 0) 


$(4,0) 
'I'(4,0; -4) 
*(4, 0: -2) 
*(4,0:2) 


ITi.Ti.T) 


5 (1, 1) (i i) 


$(5, 1) 
*(5, 1; -4) 
*(5, 1: -2) 


iTi,Ti,n) 


6 (0, 0) (0, 0) 


$(6,0) 
*(6,0; -6) 
^(e.O; -4) 
^(e.O; -2) 



Table I: Electronic eigenstates of the Ceo"" molecular ion, 
and the corresponding slave-boson amplitudes: the latter are 
selected in order to preserve, at the saddle-point level, the 
rotational invariance of the local Hamiltonian. 



The last column of the Table contains the independent 
slave-boson amplitudes ^{n,£) and '^{n,£; AN) associ- 
ated, according to Eq. ([75]) . to all the components of 
a given degenerate- multiplet: the total number of ampli- 
tudes required by our model is therefore 35, if we consider 
the full rotationally-invariant solution, while it reduces to 
13 if we force the system into the normal state, setting 
^'{n,e; AN) = 0. 



Lagrange multipliers 

The Lagrange multipliers A'^, A^p and 11^^, intro- 
duced in Sees. Ill Dl and III El to enforce the constraint 
equations (|12m4p . form, together with the slave-boson 
amplitudes, the set of variables on which the free-energy 
functional (|53p is defined. However, similarly to what we 
established in the case of the slave-boson amplitudes, we 
must note that the symmetries of our problem greatly 
reduce the number of independent Lagrange multipliers 
required for the solution of the model, and in the follow- 
ing we will identify the form of such variables. 

Denoting with a = (rn, a) both the orbital and spin 



indices, the constraints (|12lll4p read, at the saddle-point 
level. 



E^(r) = 1, 



qN 



Q 



af3 



{flu) 

k 



(91) 



(92) 



(93) 



where P(r) is the probability distribution defined in 
Eq. (1751) . Gl/(k, r) and F/(k, r) are the normal and 
anomalous quasiparticle Green's functions, and 



^ E^rn^rn'(n'|/I/;3|n), (94) 

Tnn' 

Q% ^ E^r„^r„'(n'|/I4|n) (95) 

Vnn' 

are defined as the normal and anomalous quasiparticle 
density matrices. As for -P(r), we can then make use of 
Eq. (Its]) in order to rewrite the left-hand side of Eqs. ([92]) 
and ([93]) directly in terms of the independent slave-boson 
amplitudes, obtaining the following expressions: 



qN 



Q^[<^,^5^,p (-1)"°, 
m + (T + i , [cr = ± 



(96) 
(97) 



where the spin and orbital indices of a are opposite to 
those of a. The specific choice we have made for the in- 
dependent slave-boson amplitudes thus reflects in a very 
simplified form of the quasiparticle density matrices, and 
it is not hard to recognize in this structure the symmetry 
properties which characterize our model, i.e., the band- 
degeneracy and the SU{2) ® 0{2>) rotational invariance. 
The quasiparticle energy matrix h(k)'^^ must then be ro- 
tationally invariant as well, in order to yield quasiparti- 
cle expectation values with the same structure of the QP 
density matrices: 

(/l//3)«'5a/3, (/I4)«<5„j(-1)''°. 

The kinetic part of h(k), namely Rte(k)R, is guaranteed 
to be SU{2) (g) 0(3) invariant, since it depends only on 
the degenerate band dispersion eQ,^(k) = e(k)(5Q/3 and 
on the slave-boson amplitudes (^rn[3',vE'] (through the 
i?-matrices) which have been properly selected in order 
to yield rotationally invariant solutions. On the other 
hand, the Lagrange multipliers matrices A and 11 are, 
in principle, two generic Hermitian and antisymmetric 
matrices, respectively, and we must therefore set 



A, 



Q/3 



AS. 



n, 



Q/3 



n<5„j (-1)"° 
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in order to guarantee the rotational invariance of the 
quasiparticle Hamiltonian. 

In the end, we are left with just three Lagrange mul- 
tipliers, A'^ , A, and H, to which we can eventually add 
the chemical potential /i if we decide to solve the model 
keeping the physical electron density fixed: in the latter 
case, in fact, the chemical potential plays the role of a 
Lagrange multiplier for the number equation 



nphys = J2(^lnaCma) = ^ n(r)P(r) 



(99) 



rather than being an external parameter as in the grand- 
canonical ensemble. 



Spectral weight and low-energy excitations 



As shown in Sec. IIIFl the expression for the quasipar- 
ticle spectral weight matrix Z, defined as 



as 

doj 



1 -1 



(100) 



is given, in terms of the slave-boson amplitudes, by 

Z = RRt (101) 



by the six-fold degenerate branches 



i?ip(k) = ±V[ek(|rpP + Kp)+A]^ + |Ap, (105) 



A 



\rp\^ + \rh\^ 
rln-irlUr -2Arprl 



\rpr + \rhr 

From Eq. (|105p we can then readily establish the expres- 
sion for the low-energy spectral gap. 



V'|Ap + i?4|A|-f (Ir.p + Kp)], 



(106) 



where we have assumed, for the free-electron dispersion, 
Ek G [— ^] , W being the uncorrelated bandwidth. It 
is important, however, to keep in mind that the onset of 
superconductivity in the system is signaled by the pres- 
ence of a non-zero order parameter ip*^^ rather than by 
the opening of a gap in the spectrum: these two quanti- 
ties, in fact, are directly proportional to each other only 
in the weak-coupHng regime (J/W ^ 1 and U < J), 
where we find the solution to be BCS-like, becoming in- 
stead disentangled for stronger values of either the elec- 
tron correlation U or the pairing attraction J. 



R[{(^}] being the matrix which relates the physical elec- 
tron operators to the quasiparticle ones (see Eqs. HH]) 
and ([30)1 for its definition). In our model, this relation 
reads 



rp[$,vI/]*£^ + (-l)"+-+2r,[<i>,vI/]/„ 



(102) 

and it can be easily recognized, as in the case of the 
quasiparticle density matrices, the specific structure of 
the normal (rp) and anomalous (r/i) terms, dictated by 
the symmetries of the problem. Inserting this relation 
in the definition of the quasiparticle weight matrix, we 
finally obtain 



1, 



(103) 



which states that, for rotationally invariant solutions to 
this model, all the electronic degrees of freedom are renor- 
malized by the same factor and do not mix each other due 
to the interaction terms. 

Besides the quasiparticle spectral weight, we can actu- 
ally extract, from the saddle-point values of the slave- 
boson fields and Lagrange multipliers, the entire low- 
energy spectrum of the system, i.e., its coherent single- 
particle excitations. They are defined as the frequency- 
poles of the physical electron propagator 

Dc(k,w) = RD/(k,w)R^ = R[w-h(k)]-iR^ (104) 

and, in terms of the saddle-point variables, they are given 



IV. ILLUSTRATION OF RESULTS 

As mentioned previously, the major strength of slave- 
boson approaches relies in the possibility of obtaining, 
with a relatively low computational effort, approximate 
analytical solutions which describe quite well, on a quali- 
tative footing, the effects of electronic interactions on the 
low-energy part of the spectrum. With such methods we 
can thus investigate the behavior of a given model over 
the entire range of variability of the parameters on which 
the model is defined. 

Using the slave-boson representation introduced in the 
previous section for the description of superconduct- 
ing fuUerides, we will here illustrate the behavior of 
the saddle-point solutions across the zero-temperature 
"phase diagram" , where the external parameters of the 
model are represented by the electron density riphys and 
the two ratios U/W and J/W, which measure, respec- 
tively, the strength of the Coulomb and Jahn- Teller in- 
teractions with respect to the kinetic bandwidth W t. 
We will primarily focus on the half-filled case, where it 
can be found the most interesting experimental feature of 
these systems, namely the relatively high superconduct- 
ing critical temperature in comparison to the strength of 
the pairing coupling, and then we will briefly analyze how 
the superconducting behavior extends to finite values of 
doping. The solutions are obtained using, for simplicity, 
a flat density of states, -D(e) = since the qualitative 
behavior of the system does not depend much on the 
specific form of -D(e). 



14 



— z 



/ A/J 




// 



I I I I I I I I I I I I I 

0.1 0.2 0.3 0.4 0.5 

JAV 

Figure 1: (Color online) Normal and superconducting solu- 
tions at half-filling and (7 = 0: Zn and Z are the correspond- 
ing quasiparticle weights, A and ipse are the superconducting 
gap and order parameter, (fi — iljv) represents the free-energy 
difference between the normal and the superconducting state. 



A. Half-filling (riphys = 3) 

We begin the analysis of the half-filled model illus- 
trating, in Fig. [U the normal and superconducting so- 
lutions obtained, at U — 0, for increasing values of the 
Jahn- Teller coupling J/W. Since we are turning off the 
Coulomb repulsion, this is a purely attractive model, with 
the Jahn- Teller coupling playing the role of an attractive 
local interaction acting on spin and orbital degrees of 
freedom; the physics of this system is therefore charac- 
terized by the competition between singlet formation and 
kinetic delocalization, and we find the results of this com- 
petition to be remarkably different whether we are con- 
sidering a purely normal-state solution or we are allowing 
for a superconducting order parameter. As expected for 
a purely attractive interaction, the superconducting so- 
lution is always energetically favored at finite J. 

As soon as the pairing interaction J is turned on, the 
behavior of the normal-state solution is initially char- 
acterized by a slow decrease of the quasiparticle weight 
Zn from the non-interacting value Zjq{Q) = Z{Qi) = 1, 
which is then followed by a steep descent towards zero 
for J/W > 0.2; finally, when the coupling is further in- 
creased, the metallic state turns into an insulating one, 
where all the electrons are locked in local singlets formed 
by two or four electrons"^'^, the binding energy of the sin- 
glet configuration being much more favorable than the 
kinetic energy gain associated to the electron hopping. 
This state is analogous to the paired insulator found, 
at strong coupling, in the normal solution of the at- 
tractive Hubbard model within DMFT'*'*. On the other 
hand, if we do allow for superconducting ordering, we 
find the static singlet-formation mechanism characteriz- 
ing the normal solution to be replaced by the more fa- 



vorable Cooper pair formation, in which the singlet pairs 
can still gain some kinetic energy through their propaga- 
tion: the solution, in this case, is therefore characterized 
by a finite quasiparticle weight Z over the entire range 
of the pairing interaction. We must however notice that 
the difference in the behavior of the normal solution be- 
tween the weak and strong coupling regimes (metallic vs. 
insulating) can still be traced in the behavior of the su- 
perconducting solution. In fact, for increasing values of 
J/W, we observe a crossover between a weak-coupling 
BCS-like superconductivity, where the gap A/J and the 
superconducting order parameter i/j^^ are proportional to 
each other and exponentially small in the pairing cou- 
pling (Fig. [2D, 

A/J = - ie^i A=f(J/iy) (107) 

and a strong-coupling superconductivity associated to 
Bose-Einstein Condensation (EEC) of preformed pairs, 
where both the gap and the superconducting order pa- 
rameter are saturated'^^. While in the former regime the 
formation of Cooper pairs subtracts some kinetic energy 
from the normal state in order to gain the binding en- 
ergy associated to the superconducting singlets, as evi- 
denced by the lower spectral weight Z < Zn, which cor- 
responds to more localized particles, in the large J/W 
regime, where the local singlets are already formed, the 
energy gain of the superconducting state is due to the 
finite kinetic energy of the Cooper pairs in comparison 
to the static singlets characterizing the insulating normal 
staters.. 

The most interesting aspect of the half- filling solutions, 
however, is represented by the behavior of the quasiparti- 
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Figure 2: (Color online) Superconductivity in the weak- 
coupling regime: slave-boson solution (circles and crosses) vs. 
BCS estimate. For J/W <C 1 and [7 = the superconducting 
parameters satisfy the BCS relation A/J — and, apart 

from a constant prefactor, they follow the BCS functional 
form A(J)bcs ~ We~A (\ — ^(J/Vl^) is the dimensionless 
superconducting couphng constant). 
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cle weight and of the superconducting observables (spec- 
tral gap and order parameter) as functions of the electron 
correlation t/, at different values of the Jahn- Teller cou- 
pling J. 

In Fig. [3] we plot the [/-dependence of the normal quasi- 
particle weight and of the observables characterizing 
the complete solution {Z, A and ■0^^) at J/W = 0.04, cor- 
responding to a coupling strength located in the upper 
end of the weak-coupling regime, or, in other words, just 
before the U = crossover region. The relevant feature 
to be noticed in this figure is the non-monotonic behavior 
of the superconducting parameters for increasing values 
of the electron correlation: while at small U the net ef- 
fect of the Coulomb repulsion is to rapidly destroy the 
superconducting order, as expected in a weak-coupling 
BCS superconductor (notice the small value of the su- 
perconducting order parameter at C/ = and its sudden 
disappearance as soon as U is turned on), at larger U 
values the system turns back superconducting, with an 
enhancement in the values of A and i/'sc comparison 
to [/ = 0, until it undergoes a first-order Mott transition 
at C/ = Uc, just above the corresponding Mott transition 
of the normal state. It is evident the huge enhancement 
of the superconducting amplitude i/)^^ with respect to the 
non-correlated regime. 

From a physical point of view, the reemergence of su- 
perconductivity at large U has been explained in terms 
of the "strongly correlated superconductivity" scenario 
put forward using DMFT to solve the same model^-iH 
and a related simplified model^li^. The key mechanism 
is the different effect of the correlation on the various 
interaction terms: when strong repulsion freezes charge 
fluctuations, the resulting strongly correlated quasipar- 
ticles experience a strongly reduced repulsion, while the 
attraction is essentially unscreened. As a result, the net 
effect is that J/W J/{ZW)^: when the electrons 



become more localized, the relative strength of the Jahn- 
Teller interaction grows in comparison to the renormal- 
ized hopping. The precise nature of the interaction, in- 
volving orbital and spin degrees of freedom, is crucial 
in this effect, and proves the ability of our rotationally 
invariant slave bosons to properly treat every kind of in- 
teraction. The superconducting behavior in this region 
is clearly non-BCS-like, as evidenced by the non propor- 
tionality between the gap and the order parameter: A/ J 
is indeed much larger (up to ten times) than ■f/'sc ^"^^ 
its maximum is located much closer to the Mott tran- 
sition than the order parameter's one. On the other 
hand, the pairing mechanism cannot be explained within 
a purely strong-coupling BEC-like picture, since in this 
case the pairing singlets are not already "preformed" in 
the normal phase (which is either a correlated metal or 
an S = 1/2 Mott insulator) and their fraction is much 
smaller than in standard BEC superconductivity. We are 
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Figure 3: (Color online) i7-dependence of the half-filling solu- 
tions at J/W = 0.04. The condensation energy {Q — Qm)/W 
is multiplied by a factor of 10 for visibility reasons. 



Figure 4: (Color online) [/-dependence of the half-filling so- 
lutions for fixed ratios J/U = 0.02 (top) and 0.01 (bottom). 
Circles and crosses represent the superconducting parameters 
evaluated at [/ = for the same values of J. 
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Figure 5: (Color online) i7-dependence of the half-filling solu- 
tions at intermediate-to-strong pairing couplings J/W = 0.1 
(top) and 0.2 (bottom). 



rather in the presenee of a strongly correlated supercon- 
ductor, in which a small local pairing coupling turns out 
to be enhanced, rather than suppressed, by the effects of 
a strong on-site repulsion. 

The correlation-driven enhancement of superconduc- 
tivity in the proximity of the Mott transition is even 
more evident in Fig. [31 where the solutions are evaluated 
at a fixed ratio J/U for increasing values of the correla- 
tion; together with the normal and complete solutions, 
we have plotted for comparison the (BCS-like) supercon- 
ducting parameters A/ J and V'L obtained, at U = 0, 
for the same values of J. Besides the different relation 
between A/ J and ■0°^ in the correlated case compared to 
the U — solutions, these plots emphasize how the en- 
hancement of the superconducting gap becomes stronger 
(up to three orders of magnitude in the lower panel) at 
smaller values of the pairing coupling: for J/W ^ 1, 
indeed, the value of A/ J in the proximity of the Mott 
transition turns out to be 0(1), while it is exponentially 
small in J/W in the BCS regime (see Fig. [2]). 



As already found in DMFT in a two-orbital model^, 
a completely different scenario is instead observed for 
larger values of the Jahn- Teller coupling, corresponding 
to the strong-coupling regime of the U = attractive 
model (shown in Fig. [5]). In these cases, in fact, the su- 
perconducting order parameter decreases monotonically 
with U from the large U = value, until a weakly first- 
order Mott transition (which becomes second-order when 
J/W is increased) turns the system into an insulator; a 
similar behavior characterizes the superconducting gap, 
except for an initial rise at small values of U/W in the 
case J/W = 0.1. At strong-coupling values of the pair- 
ing attraction the superconducting solution is therefore 
always energetically-favorable compared to the metallic 
one, and the electronic correlation has only the effect of 
reducing, throughout the non-insulating phase, the su- 
perconducting ordering. 

We conclude the analysis of the half-filling solutions 




Figure 6: (Color online) J-dependence of the half-filling solu- 
tions in the strongly-correlated regime: the top panel corre- 
sponds to a fixed value of the correlation, U/W = 2.45, while 
in the bottom panel the solution follows the Mott-transition 
line C/c(J). 
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showing, in Fig. [6l the non-monotonic behaviors of the 
superconducting parameters A/ J and ip*^^, as functions 
of J, in the strongly-correlated region of the phase dia- 
gram: in the top panel the value of U/W is held fixed, 
while in the bottom one it follows the Mott-transition 
line from below, U{J) = Uc{J) — SU. Combining these 
results with the ones discussed previously, we can then 
infer the existence of a second region in the J-U plane, 
beside the strong-coupling BEC-like region at J ^ U, 
in which superconductivity is found to be optimal: al- 
most surprisingly, it is located at very small values of the 
pairing attraction J and at correspondingly large values 
of the on-site electron repulsion U, just before the Mott 
localization transition line. 

The results presented in this section confirm that the 
rotationally invariant slave boson approach is able to ac- 
curately treat interactions of different kinds and, partic- 
ularly, it is not limited to charge interactions. Indeed we 
have found that the present approach is able to reproduce 
the relevant physics of a three-orbital Hubbard model for 
the fuUerenes, and, in particular, the huge enhancement 
of phonon-mediated superconductivity in the proximity 
of the Mott transition. The only qualitative aspect of the 
DMFT solution which is not found in the present study 
is the second-order (or very weakly first-order) character 
of the superconductor-insulator transition. 



B. Finite doping (iiphys = 3 — (5) 

In this section we consider the effect of doping the half- 
filled three-band model. While this situation can not be 
experimentally realized, at the moment, in fuUerides, the 
effect of doping is clearly suggestive for analogies with 
the physics of cuprates. 

The behavior as a function of doping, in the neigh- 
borhood of riphys = 3, is shown in Fig. [7] for correla- 
tion strengths respectively below and above the critical 
Mott-transition value Uc{J)', in both cases, however, we 
have U ^ J, so that they both belong to the strongly- 
correlated region of the phase diagram, where the pres- 
ence of a finite superconducting order parameter is due 
to the localization-driven enhancement of the effective 
pairing coupling. 

We find that for U <Uc the superconducting parame- 
ters decrease monotonically upon doping (eventually van- 
ishing at larger doping values), while the normal and su- 
perconducting quasiparticle weights increase from their 
finite half-filling values; as long as the superconducting 
order is present, we have Z < Zjf. On the other hand, 
for U > Uc we observe a dome-shaped behavior in the 
superconducting parameters, the gap reaching its maxi- 
mum at a very small doping value ^opt ~ 0.1, while the 
order parameter being maximum at a slightly larger opti- 
mal doping Jqp' « 0.15. In this case, both the normal and 
superconducting quasiparticle weights grow linearly with 
the doping S, but while in the overdoped region S > Sopt 
we find the standard weak-coupling behavior Z < Zjv, at 
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Figure 7: (Color online) Doping dependence of the solu- 
tions for correlation strengths respectively below (top panel, 
U/W = 2, J/W = 0.03) and above (bottom panel, U/W = 5, 
J/W = 0.02) the "■phys — 3 Mott transition. 



lower dopings we have Z > Zpj. 

The behaviors of both the normal and superconduct- 
ing solutions in the two correlation regions U ^ Uc are 
therefore remarkably different from each other; however, 
at a closer sight, we find that they can be actually ex- 
plained through the same physical mechanism, namely 
the competition between Mott-localization, which can 
eventually enhance the superconducting pairing as we 
have seen in the discussion of the half-filling solutions, 
and the delocalization tendency introduced by doping. 
In fact, when the correlation strength at half- filling is 
not large enough to completely destroy the quasiparticle 
coherence (in other words, when the quasiparticle degrees 
of freedom are nor completely frozen), we find the super- 
conducting parameters to be maximum at zero-doping. 
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where the electrons are more locahzed and the enhance- 
ment of the effective pairing coupHng is stronger. When 
the system is in the Mott-insulating phase, on the con- 
trary, there are no available quasiparticles at half-filling 
in order to develop a superconducting order parameter: 
the rcintroduction of quasiparticle coherence due to a fi- 
nite level of doping becomes then essential in order to 
recover a superconducting solution. At small dopings, 
therefore, the superconducting ordering increases, due 
to the regained coherence of quasiparticles and a still 
strong enhancement of pairing due to Mott localization; 
for larger values of doping, instead, the loosening of the 
localization-induced pairing enhancement disfavors the 
superconducting ordering, which turns to decrease as in 
the U < Uc scenario. It is interesting to note that in 
the underdoped region we have Z > Zjv, which means 
that the formation of superconducting pairs is energeti- 
cally more favorable also from the kinetic point of view, 
compared to the normal state. 

V. CONCLUSIONS 

We have generalized to superconducting solutions (al- 
lowing for the spontaneous breaking of charge symme- 
try) the rotationally invariant slave boson approach in- 
troduced by Lechermann et al^ on the basis of the work 
by Li and Wolflei. The crucial ingredient of the rotation- 
ally invariant version of slave boson methods is that the 
boson fields cannot be simply seen as probability ampli- 
tudes for the different quasiparticle states, but they are 
expressed as a non diagonal density matrix that connects 
the different quasiparticle states to the whole set of phys- 
ical states. This is easily generalized to include matrix 
elements between states with different number of parti- 
cles which allow to describe superconducting ordering. 

After a thorough description of the formalism, we ap- 
ply the method to solve a three-band model which has 
been proposed and studied to understand the properties 
of alkali-doped fuUerides^i^i^i. This model has been 
previously solved using DMFT— for integer fillings, pro- 
viding a striking realization of strongly correlated super- 
conductivity, i.e., of a situation in which strong electron- 
electron correlations favor superconductivity. A crucial 



element of the model is that the pairing attraction, which 
can be modelized as an inverted Hund's rule term, only 
involves spin and orbital degrees of freedom, which are 
not heavily affected when the charge degrees of freedom 
are frozen by the proximity to the Mott transition. This 
interorbital nature of the pairing interaction is crucial to 
give rise to the correlation-driven enhancement of pair- 
ing. In this light, this model represents an ideal test 
field for our approach, which is tailored to properly treat 
interorbital interactions that cannot be expressed in a 
density-density form. Indeed the method provides good 
results for this model, and it is actually the first mean- 
field-like approach able to reproduce the enhancement of 
superconductivity observed in DMFT. 

The good performance of the method is very encour- 
aging in view of other applications. The most chal- 
lenging direction is obviously the study of the two- 
dimensional Hubbard model, which is believed to be the 
basic model to understand high-temperature supercon- 
ductivity in the cuprates. While the full solution of 
the model on a lattice appears too cumbersome, a vi- 
able direction is the use of cluster extensions of DMFT, 
such as the cellular-DMFT^^ or the dynamical cluster 
approximation's, where the rotationally-invariant slave- 
boson method can be used as an approximate analytical 
impurity solver for the cluster Hamiltonian. This ap- 
proach has been used, for example, in Refs. andl4ll for 
normal solutions without superconducting ordering. To 
investigate the superconducting properties of the Hub- 
bard model, on the other hand, our formalism can be 
applied either to the 2x2 plaquette, where it can be 
used to better understand the outcomes of fully numeri- 
cal solutions''^, or to slightly larger clusters, such as small 
rectangles in CDMFT or the 5-site "cross" in DCA, which 
have been proposed as ideal compromises between rea- 
sonable cluster size and adequate accuracy^. 
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